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Abstract. The analysis of an INTEGRAL A02 observation of the error contours of the EGRET source 3EG J1621 + 8203 is 
presented. The only source found inside the error contours for energies between 20 and 30 keV at 5cr detection significance is 
the FR I radio galaxy NGC 6251. This supports the identification of NGC 6251 with 3EG J1621 + 8203. The observed flux is 
higher and softer than observed in the past, but consistent with a variable blazar-like spectral energy distribution. 
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1. Introduction 

One of the long standing problems of modern high-energy 
astrophysics is to understand the nature of sources emit- 
ting the highest energy y-rays. The Energetic Gamma-Ray 
Experiment Telescope (EGRET) Catalog, obtained from high- 
energy y-ray observations (£ > 100 MeV) performed between 
1991 and 1995 (Hartman et al. 1999), changed our understand- 
ing of the y-ray sky. The third version of this catalog contains 
271 point sources and about 150 of them are still unidentified. 
The search for the counterparts of these sources is particularly 
challenging since the point spread function (PSF) of EGRET 
is about 6° at 100 MeV (FWHM) and the error contours are 
large, typically extent 0.5° - 1°. During recent years, much ef- 
fort has been spent in gathering all the useful data from the 
online multiwavelength archives, searching for possible coun- 
terparts. Despite some interesting results, the quest is still open 
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(for a review see, e.g., Caraveo 2002, Grenier 2004, Mukherjee 
& Halpern 2004, Thompson 2004). 

The EGRET sources of extragalactic origin are mainly 
blazars. The correlation of EGRET error boxes with flat- 
spectrum radio sources has been revisited using deeper radio 
interferometric surveys at 8.4 GHz and improved radio spec- 
tral measurements in the northern and southern sky, down to 
6 > -40° (Sowards-Emmerd et al. 2003, 2004). These in- 
terferometric data have allowed us to discriminate more effi- 
ciently between sources with a flat spectrum core from those 
with softer lobe emission. Combined with the associations pro- 
posed by Mattox et al. (2001) using older surveys at 5 < -40°, 
a total of 128 AGN counterparts have been proposed. 60% of 
the identifications are considered very likely. Five associations 
with nearby Fanaroff-Riley type I and II radiogalaxies (FRI 
and FRII, respectively), namely 3EG J0459 + 3352 (Sowards- 
Emmerd et al. 2003), 3EG J 1324 - 4313 (Cen A, Steinle et 
al. 1998), 3EG J 1646 - 0704 (Sowards-Emmerd et al. 2004), 
3EG J1736 - 2908 (Di Cocco et al. 2004), 3EG J1735 - 1500 
(Combi et al. 2003) and 3EG J1621 + 8203 (Mukherjee et 
al. 2002) open the exciting possibility of studying misaligned 
blazars in gamma rays. While blazars can be powerful y-ray 
emitters because of relativistic beaming and amplification, the 
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Fig. 1. IBIS/ISGRI deconvolved significance map of the 
3EG J 1621 + 8203 sky region (exposure of 424 ks), in the 
20 - 30 keV energy band, with a S/N threshold set to 3cr. The 
3EG J1621 + 8203 probability contours at 95% and 99% con- 
fidence level are shown. "A" indicates the unknown excess at 
a: 3.5cr. See the text for more details. 

identification of y-ray radiogalaxies poses interesting prob- 
lems. 

According to the unified scheme of AGN by Urry & 
Padovani (1995), the transition from blazars to radio galaxies 
is based on a combination of orientation and relativistic beam- 
ing. Specifically, radio galaxies of type FRJ and BL Lac ob- 
jects are similar AGN seen with jet directions at large angles to 
the observer's view or close to it. On the other hand, FR II ra- 
dio galaxies and flat spectrum radio quasars (FSRQ), share the 
same orientation. The detection by EGRET of copious y-rays 
from blazar jets implies that FR I and FR II radio galaxies are 
also intrinsically powerful y-ray emitters, but their observa- 
tion at large aspect angle should yield a much fainter flux and 
softer spectrum. This is indeed the case for Cen A (Steinle 
et al. 1998), that has a photon spectral index in the EGRET 
energy range {E > 100 MeV) of F = 2.40 + 0.28, that ap- 
pears to be slightly softer than the average value for blazars 
(F = 2.15 + 0.04), although still consistent within the measure- 
ment errors. Cen A has a large inclination angle (~ 70°), so that 
it is difficult to explain the detection of emission in the EGRET 
energy range, but some explanations have been proposed (see 
the discussions in Sreekumar et al. 1999, and Mukherjee et al. 
2002). 

Therefore, it is very interesting to search for other radio 
galaxies that may be counterparts of EGRET sources. Their 
identification with EGRET sources would have a large impact 
yielding, for instance, valuable constraints on the collimation 
of y-ray emission from the jet and the energy of the jet parti- 
cles. In addition, it would have a significant impact on popula- 
tion studies since the density of radio galaxies is far higher than 
that of blazars, and a significant reassessment of the origin of 
the difliise extragalactic y-ray background would be required. 

Among the different cases available, we focus our attention 
on 3EG J1621 H-8203, that Mukherjee et al. (2002) have associ- 
ated with the FR I radio galaxy NGC 6251 (z = 0.02488). They 
analyzed data from X-ray satellites (ROSAT, ASCA) and by 



cross-correlating with the NRAO VLA Sky Survey (NVSS), 
they found that, among the several sources found in the EGRET 
error contours, NGC 625 1 is the most promising candidate to 
be a high-energy y-ray emitter. However, it should be noted 
that inside the EGRET error contour there are several other X- 
ray and radio sources, some of them unidentified. Therefore, 
NGC 625 1 was considered the best candidate among the known 
sources. 

In order to resolve this doubt, we requested for a long obser- 
vation with the IBIS telescope (Ubertini et al. 2003) on board 
the INTEGRAL satellite (Winkler et al. 2003). IBIS has a large 
field of view of 19° x 19° at half response, with a high an- 
gular resolution of 12' sampled in 5' pixels in the low energy 
(0.015 - 1 MeV) detector ISGRI (Lebrun et al. 2003) and in 10' 
pixels in the high energy (0.175-10 MeV) layer PICsIT (Di 
Cocco et al. 2003). The point source location accuracy (PSLA) 
of ISGRI can be down to 1' for a 30cr detection (Gros et al. 
2003). These characteristics make IBIS a valuable instrument 
to search for hard X-ray (E > 20 keV) counterparts in the large 
error contours of the EGRET sources. 

Here we report the analysis of the region containing the 
3EG J 1621 + 8203 error contours and the hard X-ray emission 
from NGC 6251. 

Throughout the paper we adopted Ho - 70 km s ' Mpc 

2. INTEGRAL data analysis 

INTEGRAL observed 3EG J1621 -H 8203 from June 17, 2004, 
18''28™ to June 23, 2004, 08M3™ UTC with a 5 x 5 dither pat- 
tern. The region of interest was about 4.6° off axis, because 
the observation of the EGRET source (Prop ID 0220018, PI 
Foschini) was amalgamated with the observation of the galaxy 
cluster Abell 2256 (Prop ID 0220020, PI Fusco Femiano), 
which was the on-axis target. 

The INTEGRAL data analysis described in the present work 
was done with the most recent version of the Offline Standard 
Analysis (OSA 4 . l') whose algorithms for IBIS are described 
in Goldwurm et al. (2003). 

The spectrometer SPI (Vedrenne et al. 2003) was perform- 
ing an annealing and therefore no data are available. In addi- 
tion, given the off-axis position, the source was not visible in 
the field of views of the X-ray monitor JEM-X (Lund et al. 
2003) and the optical monitor OMC (Mas-Hesse et al. 2003), 
both also aboard INTEGRAL. 

The analysis was performed on all the available data for an 
effective exposure toward 3EG J1621 + 8203 of 424 ks. The 
only known source detected (signal-to-noise ratio S /N - 5cr) 
inside the EGRET error contours is NGC 625 1 , found at coor- 
dinates (J2000) a = 16*'33™10^ and 6 = H-82°35'46" with an 
uncertainty of 5' at the 90% confidence level (Fig.0. The count 
rate in the 20 - 30 keV energy band is 0.1 1 + 0.02 (Icr error) 
corresponding to a flux^ of (6.7 +2.1)x 10"'^ erg cm"^ s"' (90% 

' Available through the INTEGRAL Science 
Data Centre (Courvoisier et al. 2003) at 
http : / /isdc . unige . ch/index . cgi?So£t+download 

- The conversion of count rates in physical units has been done by 
normalizing to the count rates of the Crab Nebula, obtained in sev- 
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Table 1. Fluxes in the 20 - 30 keV energy band obtained in 
different epochs. In the case of ASCA and XMM-Newton, the 
value is extrapolated from the best fit model (see the notes for 
details). Columns: (1) Satellite name; (2) Date of the observa- 
tion [DD-MM-YYYY]; (3) Photon index; (4) Flux [lO''^ erg 
cm"^ s ']. The uncertainties in the parameters are at the 90% 
confidence level. 



Detector 




Date 




r 


F 


(1) 




(2) 




(3) 


(4) 




28 


-10- 


1994 


1.8 ±0.2 


0.50 ± 0.06 


BeppoSAX'" 


19 


-07- 


2001 


1.79 ± 0.06 


1.73 ± 0.07 


XMM-Newton" 
INTEGRAL'^ 


26 
17 


-03 - 
-06- 


2002 
2004 


1 q 1+0.08 

^•^^-0.05 

2.1= 


1 iq+o.08 

^■^^-0.06 

6.7 ±2.1 



" Sambruna et al. (1999). Flux extrapolated from the fit in the 0.6 - 
10 keV energy band. 

Guainazzi et al. (2003), Chiaberge et al. (2003). Flux extracted 
from the joint fit LEGS, MEGS, PDS, in the 0. 1 - 200 keV energy 
band. 

" Gliozzi et al. (2004). Flux extrapolated from the fit in the 0.4 - 
10 keV energy band. 
Present work. 

= The conversion from rate to flux was done by assuming a Crab- 
like spectrum. 



confidence level). There is no detection in higher energy bands, 
with an upper limit (3cr) of 0.05 counts s"' in the 30 - 40 keV 
energy band, corresponding to a flux of 4 x 10"'^ erg cm"^ s"' 
assuming a Crab-like spectrum. 

Although the present statistics prevent a detailed spec- 
tral analysis, we performed some simulations with xspec 
(v. 11.3.1) and the IBIS/ISGRI RMF/ARF matrices available 
in the OSA 4 . 1 to understand the flux differences between 
INTEGRAL and BeppoSAX observations. The flux and photon 
index values observed by BeppoSAX corresponds to a count 
rate of about 0.03 c/s in the IBIS/ISGRI detector (20 - 30 keV 
energy range), to be compared with the value of 0. 11 + 0.02 
from the present observation. This confirms that the flux varia- 
tions are real. 

Inside the EGRET probability contour there is only one 
other excess (marked with "A" in Fig. ^ with S /N ~ 3.5cr, 
but it does not correspond to any known source or even the 
unidentified X-ray sources found by Mukherjee et al. (2002). 
The closest source is 1 RXS J 1 63 1 2 1 .2 H- 8 1 264 1 , 8' distant from 
the excess A. The low S /N is not sufficient to claim a signif- 
icant source detection, but it is rather likely to be a residual 
artifact of the deconvolution process. 

The other detector layer of IBIS, PICsIT, had an overall 
exposure of 302 ks, but no detection was found, with an upper 
Hmit (3cr) of 1.2 X 10"'' erg cm^^ s"', in the 672 - 1036 keV 
band in the 3EG J1621 -H 8203 region. 
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Fig. 2. Spectral energy distribution of NGC 6251 
from Chiaberge et al. (2003) with the addition of the 
INTEGRAL/lSGRl data (horizontal bars are enlarged 
for clarity). The dashed line indicates the model fit to 
BeppoSAX data; the continuous line represents the model fit to 
INTEGRAL data. The new model parameters are as follows: 
source size R = 3 x lO'^ cm, injected electron luminosity 
Linj = 6.7 X lO'*'* erg/s, slope of the injection p = 1.9, 
Tmin = 10^ Tmax = 4 X 10^, magnetic field B = 0.04 G, 
escaping time fgsc - l5R/c, beaming factor S - 3.8. See 
Chiaberge et al. (2003) for more details on the SSC model 
parameters. 

3. Discussion 

The identification of NGC 6251 as a counterpart of 
3EG J 1621 + 8203 is a complex puzzle, and with the present 
work we add some new pieces. To date Mukherjee et al. (2002) 
studied the ROSAT (0.1 -2.4keV) and ASCA (0.6-lOkeV) ob- 
servations, and the radio survey NVSS at 20 cm. They also in- 
vestigated how the y-ray emission from a relativistic jet could 
decrease as a function of the observer's viewing angle and 
found that by assuming an angle of = 45° for the jet axis 
of NGC 625 1 , the flux emitted in the EGRET energy range is 
still detectable. 

A complement, at radio wavelengths, to the work of 
Mukherjee et al. (2002) was recently published by Sowards- 
Emmerd et al. (2003), with a survey at 8.4 GHz, and they 
found that the only reliable counterpart of 3EG J 1621 + 8203 
is NGC 6251. 

However, in addition to the spatial coincidence, it is neces- 
sary that the source is physically able to generate very high- 
energy y-rays. In this respect, another piece in the puzzle 
was given by Guainazzi et al (2003) and Chiaberge et al. 
(2003): they found that the spectral energy distribution (SED) 
of NGC 625 1 has two peaks, very similar to the typical SED of 
blazars (Fossati et al. 1998). Chiaberge et al. (2003) found also 
an upper limit to the jet direction of 6 < 18° from the fit of the 
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SED with the synchrotron self-Compton (SSC) model (Fig.|3- 
On the other hand, Jones & Wehrle (2002) found 6 < 47° with 
the Very Long BaseHne Interferometer (VLBI) measurement 
of the jet-counterjet ratio. Moreover, Chiaberge et al. (2003) 
evaluated also the beaming factor 6 - 3.2, that is lower than 
for blazars (8 < (5 < 23), but greater than the case of Cen A 
(5= 1.2). 

The SSC model establishes a link between the radio and 
the hard X-ray emission, but before the present work, the 
only observation at energies greater than 10 keV was per- 
formed by the BeppoSAX PDS. However, there could be some 
doubt about the flux evaluation with the PDS, since it is a 
collimated detector with a field of view (FOV) of 1.3° and 
there is an active Seyfert nucleus at a distance of 27' from 
NGC 6251 (see the source n. 17 in the Fig. 3 of Mukherjee 
et al. 2002). Thanks to the unprecedented angular resolution 
of IBIS/ISGRI, it is possible to distinguish between the dif- 
ferent contributions. The Seyfert nucleus is not detected in the 
present INTEGRAL observation and we can exclude any con- 
tamination. Moreover, IBIS/ISGRI shows that no other source 
or excess is present in the error contours of EGRET down to 
a flux of X! 5.4 X 10"'^ erg cm"^ s"' (4cr, 20 - 30 keV energy 
range). 

The SED of NGC 6251, updated with the data from the 
present observation, is shown in Fig. |2l The model SED has 
been derived using the numerical code described in Chiaberge 
& Ghisellini (1999). The data at different wavelengths are not 
obtained simultaneously, therefore, it is not possible to draw 
firm conclusions. However, with the INTEGRAL data it is pos- 
sible to fit the data with the SSC, showing the source in a state 
with an overall higher flux. The parameters of the fit are similar 
to those in Chiaberge et al. (2003): in the present case the size 
of the source is slighlty smaller by about a factor of less than 2 
and the beaming factor increases from 5 = 3.2 to 3.8. Although 
the non-availability of simultaneous data prevents us from per- 
forming a detailed spectral modeling, we speculate that the 
"high state" observed during the INTEGRAL pointings might 
be produced by a smaller region, located even closer to the jet 
base. 

Another key question is still unexplained: the X-ray vari- 
abiUty. It is indeed known that blazars display strong variability 
on dififerent time scales, from hours to days (see e.g. Wagner & 
Witzel 1995, UMch et al. 1997). Table^shows the flux values 
obtained from recent observations with X-ray satellites. There 
are clear indication for both flux and photon index variabilities. 

EGRET observed flux variations of 3EG J1621 -i- 8203 (see 
Nolan et al. 2003). The source is not persistent, i.e. it has been 
detected in individual EGRET observations, but not in the cu- 
mulated 4-yr data. 

4. Conclusions 

The findings of this INTEGRAL observation can be summa- 
rized as follows: NGC 6251 is the only hard X-ray (20 - 30 
keV) source inside the EGRET error contours, down to a flux 
of 5.4 X 10"'^ erg cm"- s"' (4cr). The flux of the present ob- 
servation is higher than was measured by BeppoSAX (and than 
the extrapolations into the hard X-ray band of the ASCA and 



XMM-Newton spectra), thus confirming that NGC 625 1 is vari- 
able. The spectral energy distribution can be still modeled by 
the SSC model, although with slightly different parameters, but 
in agreement with a blazar-like behaviour. 

These are two more pieces added to the puzzle of the iden- 
tification of the EGRET source 3EG J1621 + 8203, but the final 
piece will be provided by the GLAST satellite, NASA's next 
very high-energy y-ray mission scheduled for launch in 2007, 
whose characteristics will assure a firm decision on the associ- 
ation of the high-energy y-ray source 3EG J 1621 + 8203 with 
NGC 6251. 
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